function [LoadRaw, RenProd, PPRegion, PPRegion2, PPInt, PPCap, PPInf, PPScale, ...
 ProdMedian,DamagesCoef1,DamagesCoef2,DamagesCoefPol1,DamagesCoefPol2]=read_elec_data()

global nnerc nhours npp npol pol_types scc

LoadRaw=zeros(nhours, nnerc);
RenProd=zeros(nhours, nnerc);


title='../Data/CleanedData/PPCap.csv';
PPCap=readmatrix(title);

npp=size(PPCap,1);




title='../Data/CleanedData/PPInt.csv';
PPInt=readmatrix(title);

title='../Data/CleanedData/PPScale.csv';
PPScale=readmatrix(title);

PPRegion=zeros(npp, nnerc);

PPRegion2=zeros(npp,nnerc);



PPInf=zeros(npp,nnerc);

%DEBUG=1
for r_i=1:nnerc
    title=(sprintf('../Data/CleanedData/LoadRaw_r%d.csv',r_i));    
    LoadRaw(:, r_i)=readmatrix(title);
    
    title=(sprintf('../Data/CleanedData/RenProd_r%d.csv',r_i));    
    RenProd(:, r_i)=readmatrix(title);   
    
    title=(sprintf('../Data/CleanedData/PPRegion_r%d.csv',r_i));    
%    title=(sprintf('../Data/CleanedData/PPRegion-own-region_r%d.csv',r_i));    
    PPRegion(:, r_i)=readmatrix(title);   
    
    title=(sprintf('../Data/CleanedData/PPRegion2_r%d.csv',r_i));  
 %   title=(sprintf('../Data/CleanedData/PPRegion2-own-region_r%d.csv',r_i));  
    PPRegion2(:, r_i)=readmatrix(title);     
    
    title=(sprintf('../Data/CleanedData/PPInf_r%d.csv',r_i));  
%    title=(sprintf('../Data/CleanedData/PPInf-own-region_r%d.csv',r_i));  
    PPInf(:, r_i)=readmatrix(title);         
    
end



title='../Data/CleanedData/ProdMedian.csv';
ProdMedian=readmatrix(title);


title=(sprintf('../Data/CleanedData/DamagesCoef1_scc%d.csv',scc));  
%title='../Data/CleanedData/DamagesCoef1.csv';
DamagesCoef1=readmatrix(title);

% put into 1000s of dollars
DamagesCoef1=DamagesCoef1/1000;

title=(sprintf('../Data/CleanedData/DamagesCoef2_scc%d.csv',scc)); 
DamagesCoef2=readmatrix(title);

DamagesCoef2=DamagesCoef2/1000;

DamagesCoefPol1=zeros(npp, npol);
DamagesCoefPol2=zeros(npp, npol);
for p_i=1:npol
    pol_type=pol_types(p_i,1);
    
    title=(sprintf('../Data/CleanedData/DamagesCoef1_%s_scc%d.csv',pol_type,scc));  
    
    DamagesCoefPol1(:, p_i)=readmatrix(title);
    
    title=(sprintf('../Data/CleanedData/DamagesCoef2_%s_scc%d.csv',pol_type, scc));  
    
    DamagesCoefPol2(:, p_i)=readmatrix(title);    
end
    
DamagesCoefPol1=DamagesCoefPol1/1000;
DamagesCoefPol2=DamagesCoefPol2/1000;

%DamagesCoef1=sum(DamagesCoefPol1,2); % do this because we corrected for NAs in data clean
%DamagesCoef2=sum(DamagesCoefPol2,2);




